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We prove the existence of a new Hamiltonian that can be used to study strongly correlated matter, 
which consists of the total energy at temperature equals zero (Eq) and the ionization energy (^) as 
eigenvalues. We show that the existence of this total energy eigenvalue, Eo does not violate the 
QQ , Coulombian atomic system. Since there is no equivalent known Hamilton operator that corresponds 

. quantitatively to ^, we employ the screened Coulomb potential operator, which is a function of this 

' ionization energy to analytically calculate the screening parameter (a) of a neutral Helium atom in 

fSJ I the ground state. In addition, we also show that the energy level splitting due to spin-orbit coupling 

is inversely proportional to ^ eigenvalue. 
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I ■ I. INTRODUCTION 

bJO' Finding even an approximate but an accurate solution to a Coulombian many-body problem (many-electron atoms 

^ , and solids) is no doubt, one of the most important problems in physics P, The Rayleigh-Ritz variational (RRV) 

O ■ method Q is a common and a powerful method regardless of the type of potential used (neglecting e-e interaction) 

'c/) ! to solve a particular Hamiltonian [1, Q. In the presence of e-e interaction, one can still use the RRV method, for 

J>~>' example, in the presence of screened Coulomb potential (Yukawa-type) to calculate the total energy eigenvalue 0, [l] ■ 

[ The Yukawa-type potential @, V ~ — (1/r) exp[— ar] is one of the well studied potential and the focus is usually 

, on how to enumerate a for a given bound state O, [HI, [13, [H, H E El, [13, [H, E S [HI, [H- In this work 
however, we will define the screening parameter, a as a function that depends on the ionization potential [23 |. in 

0^ . which, we will treat cr as a many-body parameter in real atoms. On the other hand, there are also methods employed 

' to tackle screening effect namely, renormalization-group 2^ , 1/iV expansion [2^, [23 , screened Coulomb potential 



X 



m 



I the momentum representation [2g| and, Dirac and Klein-Gordon equations with relativistic effects [H, [10, [3l|. These 

■ methods, though very useful and accurate, do not allow one to associate the different isolated atomic energy levels 
\ for all the ions that exist in solids to changes on the electronic excitation probability in solids. This association is 
. important to predict electronic properties of solids with different ions. 

[ — ' As such, we propose here a new many-body Hamiltonian to study just that, using the Yukawa-type potential. 
. The eigenvalue of this Hamiltonian is made up from the total energy at temperature equals zero (-Eo) and ionization 
' energy (^). The total energy from this Hamiltonian has been verified in experiments where the isolated atomic energy 
O I levels have been used to predict the changes of the electronic excitation probability in solids (Ref. [l^l and references 
' therein). Hence, in this paper we do not discuss the application of this Hamiltonian in solids. As such, our motivation 
J>~>] here is to give proofs-of-existence that the right-hand-side (eigenvalue) and the left-hand-side (Hamilton operator) 
1-^ . of this new many-body Hamiltonian are physically equivalent. We apply this Hamiltonian to Hydrogen and Helium 
atoms. In order to do so, we will need to find an operator that corresponds to ^ eigenvalue and prove that it yields 
the same conclusion as ^ itself. We first establish the validity of this eigenvalue for the ID systems and Hydrogen 
atom. Secondly, wc will show that the screened Coulomb potential (V^c) is the operator that corresponds to ^, and 
^ , the screened Coulomb potential in atomic He is also proportional to the ionization energy. We will also derive an 

■ expression that associates the spin-orbit coupling and the ionization energy, through the screener, a. 

II. MANY-BODY HAMILTONIAN AS A FUNCTION OF IONIZATION ENERGY 

The 3-dimensional Schrodinger equation for Fermions with mass, m moving in the presence of potential, ^(r) is 
given by (after making use of the linear momentum operator, p = —ifiV^) [3^ . 
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= iE + V{r))ip. 



E denotes the total energy. 
We define, 



± C := Ski, 



V{v) 



(1) 



(2) 



such that ±^ is the energy needed for a particle to overcome the bound state and the potential that surrounds it. 
i?kin and Eq denote the total energy at V{t) = and the energy at T = 0, respectively, i.e., E'kin = kinetic energy. 
In physical terms, ^ is defined as the ionization energy. That is, ^ is the energy needed to excite a particular electron 
to a finite r, not necessarily r ^ oo. 

proof: At temperature, T = and V{v) ^ 0, H(f = -^V^^s = Eoip. Hence, from Eq. (P) with V{r) = and T = 
0; the total energy can be written as 



kin 



En 



(3) 



Remark 1: Therefore, for an electron to occupy a higher energy state, N from the initial state, AI is more probable 
than from the initial state, L if the condition, [En{> 0) — Em{> 0)] < [En{> 0) — £'l(> 0)] at certain T is satisfied. 
As for a hole to occupy a lower state M from the initial state TV is more probable than to occupy a lower state L if 
the condition, \Em{< 0) - En{< 0)| < \El{< 0) - En{< 0)| at certain T is satisfied. 

On the other hand, using the above stated new definition (Eq. ^) and the condition, T = 0, we can rewrite the 
total energy as 



E — i?kin — Eo±^. 



(4) 



Therefore, Remark 1 can be rewritten as 

Remark 2: For an electron to occupy a higher energy state, N from the initial state, M is more probable than from 
the initial state. L if the condition, ^(M — > TV) < ^(L — > N) at certain T is satisfied. As for a hole to occupy a lower 
state M from the initial state TV is more probable than to occupy a lower state L if the condition, ^(TV M) < 
^(TV — + L) at certain T is satisfied. 

For electron-like excitations. Remark 2 implies that 



■N)- 



^ [ENi> 0) - 



- L^.vv- -J -Em{> 0)], 
--[En{>Q)-El 



^(>0)], 



while 



e(TV ^ M) = \Em{< 0) - £;Ar(< 0)1, 
i{N ^ L) = \El{<^)- En{<0)1 

for hole-like excitations. In other words, if we let the energy function, E = EqzLS^ 
the many-electron atomic Hamiltonian as 



£^kin + ^(r), then we can write 



Hip = [Eo ± e)(^. 



(5) 



III. INFINITE SQUARE- WELL POTENTIAL 

Here, the hamiltonian given in Eq. ([5]) is applied to a free-electron system, by considering a free-particle of mass m 
moving in 1-dimension of an infinite square well (width = a). Therefore, Eq. ([5]) can be solved generally, to give 



Csin 



1/2 



(6) 
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After normalization, jCp sm'^[{^{Eo±^)y/'^]dx = |Cpa/2 = 1, one obtains C = Finally, the normalized 

wave function. 
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Applying the boundary conditions for free-electrons, V{x) = 0; ^ = ii^kin ~ ^^o + [Vi^) = 0] = i?kin — Eq; E = Eo±^, 
(p{0) — and 93(a) = require {^{Eq zb^))^^ = rnr/a. Note here that the condition, V{x) — that leads to the 
free-electron concept also implies the square well potential equals zero anywhere between and a (0 < x < a), and 
this will stay true for as long as a ^ 2re where denotes an electron's radius. Eventually, one arrives at 





mr 




X 




a 
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Obviously, one can also obtain the exact form of Eq. ([5|) from the ID time-independent Schrodinger equation with 
V{x) = and fc„ = mr/a (Ref. [s^). Therefore, we can conclude that our Hamiltonian is exactly the same as the 
usual Hamiltonian. However, in the presence of potential energy, our Hamiltonian will provide additional information 
on the energetics of Fermions in different atoms. This new information will tell us how one can use the atomic energy 
level difference to predict the electronic excitation probability in solids via Fcrmi-Dirac statistics (2^ . 



IV. ID DIRAC-DELTA POTENTIAL 



As a matter of fact, there can be many general solutions for Eq. ([5]) with V{x) ^ and these solutions can be derived 
in such a way that they can be compared, term by term with known wave functions. For example. If V{x) = —aS{x), 
then we need a solution in the form of ip{x) = Cexp[— iax] and the associated wave function can be derived as 

^Q^^-* = ~iaC exp[~iax], ^ q^^'^ = —a'^Lp{x). we can rewrite Eq. ([5]) to get 



d'^Lp 2m 
dx^ 



(9) 



Using Eq. dH), we find a = [^{Eq ± ^)] Therefore, Lp{x) = Cexp [ - i{^{Eo ± . Normalizing Lp gives 



l/2„ 



i=c::Hxrdx 



Ym- Hence, 



(p{x) 



;i/2 



X exp 



2m , 

-j^iEo±C) 



1/4 



2m 



(Eo ± 



1/2 -I 



(10) 



Term by term comparison between Eq. (jlOp and [3; 



, , ma 
<^(.t) = — exp 



gives the bound state energy, £'0 ± ^ = 
^))i/2a; = -ma\x\/h\ 



72^2 gj^j^gj. by equating \i'^{^iEo ± 0)] 



1/4 = ma/ft or -i{^{Eo ± 



V. HYDROGEN ATOM 



In this section, Eq. ([5|) is applied to a Hydrogen atom with Coulomb potential and subsequently, exact results of 
its energy levels are derived. The radial equation of a Hydrogen atom is given by [s^l 



4 



2m dr'^ 



+ 



u ~ Eu. 



(11) 



Equation can be rewritten using = i?o ± and can be readily solved [33l | to obtain the principal quantum 
number, n 



2n 



27Teoh\/^iEo±0 



2m 



o2 1 



(12) 
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ub denotes the Bohr radius and eq is the permittivity of space, m and H are the particle's mass and the Planck 
constant, respectively. Equation can be rewritten so as to obtain the energy levels of a Hydrogen atom 



{Eo ± On = 



2h?\4TTeo) 



= En- 



(13) 



We propose that {Eq ± On can be rewritten in terms of the standard En where, —{Eq ± On = ~En- Proof: 

-(-Bo ± 0™+i =-[En± {E,n - S™+i)] ^ Eq. @ 
= - [E,n ± (-Em - E,n+i)] ^ Eq = E.„, = grouud state 

= — [Em — Em + Em.+ l] ^ Em. < Em.+ l ± ^ — 
= -~Em. + Em — Em+1 = —Em+1 

Alternatively, one can also show that 



(14) 



-(-Bo ± Orn+l =-[Ea± {Em+1 - Em)] ^ Eq. © 
= - [Em. ± {E.m+1 - Em.)] <^ Eq = E.^ = grouud state 

= — [Em. + E.m+i — Em] ^ Em.+ l > Em .'. ± ^ + 
= ^Em. + Em — B'rn+l = ^Em^i 

Equation or (fT5|) can be used to calculate {Eq ± On- Example: 



(15) 



-{Eq ± On=l = — [Eo ± {En=l — En=l)] = — 

-{Eq ± On=2 = — [En=l — {En=l — -^,1=2)] = — -E„=2 

-(£'0 ± 0n=3 = — [En=l — {En=l — En=3)] = — -En=3 



(16) 



Consequently, energy levels for a Hydrogen atom can be written exactly, either as {Eq ± On or i?„. Alternatively, 
we can also assume a solution for the Hydrogen atom in the form of (pn=i,i=o,m=o{f,b) — Cexp[— &r^] comparable 
with [33I the standard hydrogcnic radial wavefunction, 



(pioo{r,0,<p) 



1 



-r/aE 



(17) 



Therefore, we obtain hitp{r) = InC — 6rlne, 



1 dip(r) 
ip{r) dr 



-b and 



\)!£ill_ Equation ([9]) is used to obtain 



h^'-^2m[EQ±i]. 



(18) 
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Normalization requires, 1 = \ip{r)\'^(fir = AirC^ ^^^yi , consequently, ipioo(r,b) ~ y ^ ^ l/'^-B- The 

expectation value of the momentum can be calculated as, (p) = (p{r,b){~ih{d/dr))ip{r,b)(Pr — h/ias- The 
momentum can also be written as, p = h?k^ ^ ^j2m{Ei) ± hence. 



ft / 1 

as 



i y 2m(i;o ± 

Note that Eq. (fT9|l can also be obtained from Eq. (fTSjl . Equation (fT9|) is the ground state energy for atomic 
Hydrogen, as it should be and it justifies the applicability of Eq. Now, for free many-electron atoms, it is 

neccessary to invoke Eqs. ([T^ . ([TS]) and which can be easily extended exactly to free many-electron atoms. For 
example, the eigenvalue for the atomic Mn with electronic configuration, [Ar] 3d^ 4s^ can be written as (after exciting 
one of the 4s^ electron to 4p) 



-{Eo ± e)4pi = -[Eo± {Eipi - Ei,2)] ^ Eq. © 

= - [E4s2 ± (i?4pi - i?4s2 )] ^ Eq ~ i?452 = initial state 

= — [i?4s2 -|- £"4^1 — £"432] ^ E^pi > £4^2 .'. zh — > + 

= -E^pi (20) 

Using Eq. ([2]), one can also write, ^ = i^kin ~ Eq + V{x) = £'4pi — £452, which suggests that we can describe the 
excited electron's properties from the energy level difference, before and after the excitation. All the results presented 
here, thus far are exact and straight forward. 



VI. SCREENED COULOMB POTENTIAL 



In the previous sections, we only worked on the right-hand-side (RHS) of the Hamiltonian, Eq. ([5]), which is the total 
energy eigenvalue and we did not touch the Hamilton operator (LHS). In the subsequent sections however, we will need 
to solve the Hamilton operator in order to evaluate the influence of the ionization energy on the operator side (LHS), 
whether they tell us the same story. Of course, we cannot expect to find an isolated or a stand-alone operator that could 
correspond quantitatively to the ionization energy eigenvalue on the RHS. This is true for all atoms with potential 
energy attached, including the electron-electron (e-e) interaction term that gives rise to screening. Mathematically, 
we say that this one(LHS)-to-one(RHS) (Kic(r) to ^) correspondence does not exist because ^ has been defined as a 
function of both kinetic and potential energies (see Eq. However, we can repackage the Coulombian e-e repulsion 
potentials as the screened Coulomb potential (Yukawa-type) for atoms other than Hydrogen. Therefore, we can use 
the screened Coulomb potential to evaluate how the ionization energy influences the screened Coulomb potential in 
atomic He. Before we move on, let us just accentuate an important point here. Indeed, it is true that common 
techniques of many-body theory require handling of the Hamilton operator first, and then calculate the corresponding 
eigenvalue, {E, the total energy). In this paper however, we are doing just the opposite, meaning, we have redefined 
the total energy eigenvalue, as given in Eq. and now we are in the midst of verifying its consequences on the 
Hamilton operator. Of course, this verification is indirect simply because we do not have a stand-alone operator 
that quantitatively corresponds to the eigenvalue, ^. As such, we will make use of the screened Coulomb potential 
operator, which is given by Q 



47reor Aireor 

Here, /x is the screener's constant of proportionality, while A = (127reo/e^)aB- AH we have to do now is to show that 
for atomic He, {V)^^ is inversely proportional to a. Secondly, we will also need to show that (V)^^ is proportional to the 
ionization energy, by comparing with other two-electron system. For example, when we compare atomic He with other 

two-electron ions, namely, Li+ and Be^+, then we expect, (V)^^ < (V)^^ < (V)^^ , because ^hc(54.4 — 24.6 = 29.8 
eV) < CLi+ (122.5 - 73.6 = 46.8 eV) < ^bc^+ (217.7 - 153.9 = 63.8 eV). 
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VII. HELIUM ATOM 
A. Many-body Hamiltonian and screened Coulomb potential 

Here, we will first find the expectation value for the screened Coulomb potential for atomic He. The He wavcfunction 
is the product of two hydrogenic wavefunctions (recall that we are not solving the Hamiltonian for the He atom, rather 
we are only solving the screened Coulomb potential operator) [s^ . 



<^o(ri,r2) = (y5ioo(ri)<Pioo(r2) 



-Z(ri+r2)/aE 



(22) 



Subsequently, the mutual screening between the electrons, ri and Y2 implies that both electrons have identical 
effective charge of < 2e because they screen each other. Hence, 



o2 r 



47reo 

^2 



1 



1 



|ri - r2 



iOZ^'Ei 



47reo V2aB/ (22' + aBcr)5 [2Z + asa] 



(23) 



where, |ri — r2| = ^/rf + r-j — 2rir2 cos(6'2). From Eq. (|23p. it is clear that a is inversely proportional to the 
ionization energy, whereas the screened Coulomb potential is proportional to the ionization energy. The complete 
Hamiltonian for the Helium atom can be written as (with the screened Coulomb potential as the correction term, 
ignoring fine structure and other corrections) 



Where Ho is given by [33, l32 



H = Ho + Vsc 



Zm 47reo 



2 



2 



(24) 



(25) 



In Eq. p4p . we again simplified the Hamilton operator by strictly forcing the e-e interaction into the (V)^^ as the 
correction term. Therefore, the non-relativistic many-body Hamiltonian (again, ignoring the fine structure corrections) 
can now be written as 



H 



fi2 



ineo 2 



(26) 



while its eigenvalue is exactly, equals to Eq ± ^. Firstly, the aim of Eq. p6|) is to point out that the eigenvalue, 
^ and (V^)s(, provide the same conclusions. Secondly, one can use the ionization energy to make certain quantitative 
and qualitative predictions in solids [l^l, as well as in atoms (discussed here and in the subsequent sections). The 
total-energy expectation value for atomic He can be obtained (as given below) by using Eq. (p6|) . 



{H) = 8E, + {V)^ 



(27) 



Note here that Eq. is mathematically identical with Eq. ([5]), where i?o ± ^ = 8Ei + (V)^^ and ^ physically 
corresponds to (V),^ through the screener, a as defined earlier. Our problem now is to solve Eq. (|27[) . which is given 
by 



(28) 



7 



Equations ((28)) tells us that with increased ionization energy (^) , the magnitude of the screener (u) is reduced and 

eventually, gives rise to stronger e-e interaction, {V)^^. This conclusion is understandable because small ^ implies 
weak Coulomb force between the valence electron and the atomic core electrons, which in turn gives a strongly 
screened valence electron. As a consequence, this particular valence electron interacts less with the core electrons and 
defines the limit between the strongly correlated electrons oo) and free-electron solids ^ 0). Using the RRV 
method, one can capture the exact scenario stated above with variationally determined charge strength, Ze < 2e for 
He atom [s^,!!^- Z is the atomic number. Since the ground state energy of the neutral He and H are known [ssl. [39l|. 
79.005 eV and 13.600 eV, thus we can estimate the screener for atomic He, 



as 



-79.005 + 8£;i 



2.02 X 10^ 



(29) 



The limit, lim^^oo (^)sc ^ linv^o (^}sc ~ '^^ purely due to electron-electron interaction (stronger) with zero 
screening. This value agree exactly with the one calculated in Ref. [3^, also without screening. As a consequence, 
the ionization energy theory does not violate the true physical picture of the atomic He, both quantitatively and 
qualitatively. Now, let us calculate the screened Coulomb potential for other two-electron ions, namely, Li+ and 

. . - Li+ - Bo^+ - He 

Be2+. Using Eq. in the limit cr ^ 0, we obtain, (V)^^ = 51 eV and (V)^^ = 68 eV. Hence, we have (V)^^ (34 

eV) < (V)^^ (51 eV) < {V)^^ (68 eV), as expected. Furthermore, we can countercheck whether Eq. ([28|) is derived 
without mistakes by showing that it satisfies the Hellmann-Feynman theorem [ssl . [36| , 



d{H) 



d 

mz"^ El- 
der 





OH 






d<7 


^0^ 



1 



i2Z + aBa) 



47ren 



200Z^EiaB 
' {2Z + aBof 



Q /g-o-('-i+r2) 

da \ Iri - V2\ 



(30) 



In summary, we have justified the influence of the ionization energy in Eq. ([5]) through the screeners and the 
screened Coulomb potential, applied to the Helium atom via Eq. (|28p . Equation ([5]) is technically easier to use, for 
example, since ^ is unique for each atom and using the total energy, i?o ± C ^-s one of the restrictive condition, we 
can derive the ionization energy based Fermi-Dirac statistics, which can be applied to non free-electronic solids. The 
respective distributions for electron and hole are given by [24| 



/e(i?o,e) 



1 



1 

e[S^?'-(Bo-0]/feBT 



(31) 



where, for solids, E^p' is the Fermi level at T = 0, T and fc^ are the temperature and Boltzmann constant, 
respectively. 



B. Energy level splitting due to spin-orbit coupling 



In this section, our intention is to find a qualitative expression between the spin-orbit splitting and the ionization 
energy. The spin-orbit coupling operator after incorporating the screened Coulomb potential (Eq. pip ) is given by 

^ _ S-L Ze2 ri d 1^_ 
2w?(? 47reo \r dr r 

= DrV + Vk (32) 
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where 137 



(D) 



(33) 



c is the speed of hght, S and L arc the spin and orbital angular momentum operators, respectively. The / in 
(Z, — / — 1) is for i ^ I + \, while the — / — 1 is for j = I — \. We ignore the irrelevant part of the Hamiltonian for the 
time being (labeled D) and focus on the term in the curly bracket. Here we consider the Hydrogenic wavefunction 
again, given by [3^ 



<y92io(ri,r2) 

--{- 

327r \ aB 



T\r2e COS 6^1 cos 6/2 . 



(34) 



Recall here that we can use the appropriate Hydrogenic wavefunction, depending on the type of orbital, since our 
corrections have been incorporated into the Screened Coulomb potential (see Eq. ((27)) ). From Eq. ([32]) we have. 



^2W 



</'210 



(35) 



We solve the two potential terms separately, 



(D) • V2 



210 



(D) 



(327r)2 \aB 



¥5210 
10 



r'ie "^'■1 cos^^i (fri 



' rl 



(D)r.^ 



cos"^ 6*2 d r2 
3 



B 



(36) 



and 



(D) • ( (^210 



</5210 



\aB) 24(CTas + Z)2- 



Finally, 



{H] 



soc 



as 



l2{c7aB + ZY 24{aaB + Z) 



(37) 



Apart from (V)^^, we also find that the spin-orbit separation, -ffsoc is proportional to ^. The -ffsoc ^Iso satisfies 
the Hellmann-Feynman theorem. 



d{H) 



soc 



da 





dHsoc 




da 



<^210 



-(D) 



4(0-05 + Z)'^ \aB 



(38) 
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Using the value for screener calculated earlier (2.0 x 10^ m we find a crude estimate for the spin-orbit separation 
for He in 2p-orbital, which is given by, 

HH)soc = 3-9 X 10"^(;, -I - 1) eV 

= 0.72 meV ^ / = 1. (39) 

This calculated value is reasonable compared to the relativistic corrections [s^ for Is orbital {I — 0), which is 
0.00311 eV. Apparently, <5(i/)gQQ will be smaller even for / = 1 as a result of relativistic kinetic energy exclusion, 
because relativistic correction for / = 1 must come from both spin-orbit coupling and relativistic kinetic energy (s^ . 

VIII. APPLICATION TO MANY-ELECTRON ATOMS 

There is one important relationship, given in Eq. (|37p that we can use it to qualitatively check the energy level 
splitting in many-electron atoms. Equation (P7|) says that the magnitude of the energy level splitting {5), of a 
particular orbital is proportional to the energy level difference (^), of two different orbitals. That is, if the energy 
level difference is large, then it implies that S is also large. The S and ^ can be compared to the experimental atomic 
spectra using the following approximate formulae, 

^= yE'^^""- ""El'"^, (40) 

and 

where x and y denote the different electronic configuration (EC) and/or spectroscopic term. That is, for atomic He, 
X = ls3p:^P and y = ls3d:^D, for atomic C, x = 2s^2p^:^P and y = 2s^2p^:^D, and so on. The magnitude of ^ and S 
can be calculated for atomic Mn using Eqs. (gQl) and gj): VEf'''' = 2.319170 eV, ''Ef''' ^ 2.114214 eV, therefore ^ 
= 2.319170 - 2.114214 = 204.96 meV. Whereas, = 2.186728 eV thus, S = 2.186728 - 2.114214 = 72.51 meV. 

Here, three atomic spectra namely, He-I, C-I and Mn-I, as representatives of the periodic table are considered and are 
available in Ref. However, there are also atoms where the energy level difference is not clearly defined because of 
the energy- level overlapping, and incomplete spectroscopic data (for example, atomic Nd-I). The selected electronic 
configuration (EC), J (total angular momentum), spectroscopic term, ^, energy levels {E^) and its splittings (S) are 
listed in Table [J It is clear from Table |T] that the ionization energy concept invoked in the many-body Hamiltonian 
(see Eqs. ([5]) and ()26p ) explains the energy level splitting in many-electron atoms where, the energy level splitting, S 
is proportional to the energy level difference of a particular atom. 

IX. CONCLUSIONS 

In conclusion, we have established a new many-body Hamiltonian where the screened Coulomb potential operator is 
found to correspond directly qualitatively with the ionization energy eigenvalue. It is not possible to find an operator 
that corresponds quantitatively due to the definition of the ionization energy itself. The advantage of working with this 
Hamiltonian is that we do not need to solve the Hamiltonian in order to make qualitative predictions in many-electron 
atoms. In other words, we can use the total energy, which is the function of ionization energy in order to predict the 
changes in screening strength and spin-orbit separation for different atoms. Here, we have given proofs of existence by 
applying it to one-dimensional systems, atomic hydrogen and helium, and ionic lithium and beryllium. We found that 
the ionization energy is a mathematically and physically valid eigenvalue and its corresponding screened Coulomb 
potential operator is also found to be valid. 
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TABLE I: The energy level difference (^) and the corresponding energy level splitting (5) were calculated using the experimental 
data obtained from Ref. Eqs. ^ and (|4T|) . 
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